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Modern data is showing increasing evidence that the Universe is accelerating. So far, all attempts 
to account for the acceleration have required some fundamental dimensionless quantities to be 
extremely small. We show how a class of scalar field models (which may emerge naturally from 
superstring theory) can account for acceleration which starts in the present epoch with all the 
potential parameters 0(1) in Planck units. 



PACS Numbers : 98.80. Cq, 95.35+d 

Current evidence that the Universe is accelerating [0 , 
if confirmed, requires dramatic changes in the field of the- 
oretical cosmology. Until recently, there was strong prej- 
udice against the idea that the Universe could be acceler- 
ating. There simply is no compelling theoretical frame- 
work that could accommodate an accelerating universe. 
Since the case for an accelerating universe continues to 
build (see for example Q), attempts have been made to 
improve the theoretical situation, with some modest suc- 
cess. Still, major open questions remain. 

All attempts to account for acceleration J3| in- 
troduce a new type of matter (the "dark energy" or 
"quintessence" ) with an equation of state pq = wq pq re- 
lating pressure and energy density. Values of wq < —0.6 
today are preferred by the data [p2| and in many mod- 
els wq can vary over time, (fn this framework, wq = 
— 1,wq — gives a cosmological constant.) 

One challenge faced by quintessence models is similar 
to the old "flatness problem" which is addressed by cos- 
mic inflation. Consider Qtot = Ptot/pc where p c is the 
critical density (achieved by a perfectly fiat universe). 
The dynamics of the standard big bang drive Q t ot away 
from unity, and require extreme fine tuning of initial con- 
ditions for tttot to be as close to unity as it is today (infla- 
tion can set up the required initial conditions) . In models 
with quintessence one must consider Qq = pq/ p c which 
is observed to obey 



Mother = (Ptot ~ Pq) / Pc 



(1) 



today. The a "fine tuning" problem in quintessence mod- 
els comes from the tendency for J7q to evolve away from 
Mother- Equation ^ is achieved in these models either 
1) by fine tuning initial conditions or 2) by introduc- 
ing a small scale into the fundamental Lagrangian which 
causes pq to only start the acceleration today. This sec- 
ond category includes cosmological constant models and 
also a very interesting category of "tracker" quintessence 
models 0Jll],[l8| which achieve the right behavior inde- 
pendently of the initial conditions for the quintessence 
field. One then has to explain the small scale in the La- 
grangian, and this may indeed be possible (23) 

Here we discuss a class of quintessence models which 
behave differently. Like the "tracker" models, these mod- 



els predict acceleration independently of the initial con- 
ditions for the quintessence field. These models also 
have pQ(today) fixed by parameters in the fundamen- 
tal Lagrangian. The difference is that all the parameters 
in our quintessence potential are O(l) in Planck units. 
As with all known quintessence models, our model does 
not solve the cosmological constant problem: We do not 
have an explanation for the location of the zero-point of 
our potential. This fact limits the extent to which any 
quintessence model can claim to "naturally" explain an 
accelerating universe. Recently Steinhardt |Q has sug- 
gested that Af -theory arguments specify the zero-point of 
potentials in 3+1 dimensions. Our zero point coincides 
with the case favored by Steinhardt 's argument. 

We start by considering a homogeneous quintessence 
field <f> moving in a potential of the form 



V{cf) 



(2) 



We work in units where Mp = = H = c = 1. 

The role of spatial variations in such a field has been 
studied in H,|^2,g3,BlJ. Inhomogeneities can be ne- 
glected for our purpose, which is to study the large-scale 
evolution of the universe. We assume inflation or some 
other mechanism has produced what is effectively a flat 
Friedmann- Robertson- Walker universe and work entirely 
within that framework. Cosmological fields with this 
type of exponential potential have been studied for some 
time, and are well understood p5|-|3"o||. (A nice review 
can be found in |3l[] . ) 

Let us review some key features: A quintessence field 
with this potential will approach a "scaling" solution, in- 
dependent of initial conditions. During scaling Qq takes 
on a fixed value which depends only on A (and changes 
during the radiation-matter transition). In general, if 
the density of the dominant matter component scales as 
p oc a~" then after an initial transient the quintessence 
field obeys = , effectively "locking on" to the dom- 
inant matter component. Figure |l|(upper panel) shows 
Qq(ci) for scaling solutions in exponential potential mod- 
els, where a is the scale factor of the expanding universe 
(a(today) = 1). At the Planck epoch a ps 1CT 30 . In Fig. 
|l|( upper panel) one can see the initial transients which 
all approach the unique scaling solution determined only 
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by A. In |32] it is shown that exponential potentials are 
the only potentials that give this scaling behavior. 

Scaling models are special because the condition Qq « 
Mother is achieved naturally through the scaling behavior. 
The problem with these exponential models is that no 
choice of A can give a model that accelerates today and 
is consistent with other data. The tightest constraint 
comes from requiring that J7q not be too large during 
nucleosynthesis [[33| (at a « 1CP 10 ). The heavy curve in 
Fig [l](uppcr panel) just saturates a generous < 0.2 
bound at nucleosynthesis, and produces a sub-dominant 
J7q today. The combined effects of sub-dominance and 
scaling cause wq = in the matter era, so this solution 
is irrelevant to a Universe which is accelerating today. 
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FIG. 1. The upper panel shows £Iq (a) for different constant 
values of A. Each solution shows scaling behavior after the 
initial transient. The radiation-matter transition is evident 
at a ~ 10~ 5 . The heavy curve saturates a generous interpre- 
tation of the nucleosynthesis bound, but still does does not 
generate acceleration today. The lower panel shows the evo- 
lution of 4>{a) for the same solutions. Note how 4> varies very 
little while a and pQ vary by many orders of magnitude. To- 
day a — 1, a ~ 10~ 10 at nucleosynthesis, and a ~ 10 -30 at the 
Planck epoch. The curves correspond to A = 1.5, 2.5, 3.5, 4.5, 
and 5.5 (going from top to bottom in the upper panel). 

The lower panel of Fig. [I] shows how the value of cj> 
changes by only about an order of magnitude over the 
entire history of the universe, while the scale factor (and 
Pq) change by many orders of magnitude. This effect, 
which is due to the exponential form of the potential, 
plays a key role in what follows. The point is that mod- 
est variations to the simple exponential form can produce 
much more interesting solutions. Because <j) takes on val- 
ues throughout history that are O(l) in Planck units, the 



parameters in the modified V(4>) can an also be 0(1) and 
produce solutions relevant to current observations. 

Many theorists believe that fields with potentials of the 
form 



V{4>) = V p (0)e- 



(3) 



are predicted in the low energy limit of M-theory , where 
V p (4>) is a polynomial. As a simple example we consider 



V p {<t>) = (0 - B) a + A. 



(4) 



For a variety of values for a, A and B solutions like the 
one shown in Fig. ^| can be produced. In this solution 
J7q is well below the nucleosynthesis bound, and the uni- 
verse is accelerating today. We show ^Iq(o) (dashed), 
^matter («) (solid) and Q m diaUon(a) (dotted). The lower 
panel in figure || plots wq (a) and shows that the neces- 
sary negative values are achieved at the present epoch. 
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FIG. 2. Upper panel: A solution obtained by including a V p 
factor in the potential. We show Oq(ci) (dashed), f2 maH(!r (a) 
(solid) and Eradiation (a.) (dotted). The lower panel shows 
wq(o). The solution shows the normal radiation and mat- 
ter dominated epochs and then an accelerating quintessence 
dominated epoch at the end. We have used V p given by Eqn. 
gwith B = 34.8, a = 2, A = .01, and A = 8. Today a = 1, 
a ~ 10 -10 at nucleosynthesis, and a ~ 10~ 30 at the Planck 
epoch. 

Figure |^ illustrates how the solutions depend on the 
parameters in V p ((f>). We plot quintessence energy den- 
sity pq as a function of the scale factor a. After showing 
some initial transient behavior each solution scales with 
the other matter for an extended period before pq comes 
to dominate. The radiation-mater transition, which oc- 
curs at around a = 10~ 5 can be seen in figure a as a 
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change in the slope in the scaling domain. The constant 
parameter B in Eqn. |] has been selected from the range 
14 — 40 for these models, yet the point of <j> domination 
shifts clear across the entire history of the universe. In 
this picture, the fact that £Iq is just approaching unity to- 
day rather than 10 10 years ago is put in by hand, as is the 
case with other models of cosmic acceleration. Our mod- 
els are special because this can be accomplished while 
keeping the parameters in the potential O(l) in Planck 
units. Although we only illustrate the B dependence 
here, we have found that similar behavior holds when 
other parameters in V p (<p) are varied. 



This varies rapidly near <f> = B and for a = 2 peaks at a 
value V p = l/VA. The upper panel of Fig. || shows the 
behavior of V near <fi = B = 34.8 for the solution shown 
m Fig |. The dashed curve shows a pure exponential 
for comparison. The lower panel shows the curves V^/V p 
and A (the constant curve). Where these two curves cross 
V = 0. Because the peak value 1/yA > A, two zeros are 
produced in V' creating the bump shown in the figure. 
In our solution pq is coming to dominate near <f> = B 
because the field is getting trapped in the local minimum. 
The behavior of the scaling solution ensures that <f> gets 
stuck in the minimum rather than rolling on through 
(regardless of the initial conditions). 
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FIG. 3. Energy density (pq) vs a. The heavy shows the 
solution from Fig ^| This set of solutions shows that the 
point where acceleration (pq = const) starts is controlled by 
Planck-scale parameters in the Lagrangian. Moving from top 
to bottom on the right side the values of B (in Eqn. 4) are 
14, 22, 28, 34.8, and 40. We also varied the initial value of <f> 
between 0.3 and 10 to illustrate a range of initial transients. 



Let us examine more closely what is going on: 
derivative of V(<p) is given by 



dV 
~d4 



V 



- A V. 



The 



(5) 



In regions where V p is dominated by a single power-law 
(j) n behavior VS/V P « n/4> which, unless n is large, rapidly 
becomes << A for values of A large enough to evade the 
nucleosynthesis bound, leading to little difference from 
a simple exponential potential. However, there will be 
points were V p can show other behavior which can impact 
V. Using equation [| gives 



V 
' p 

Vn 



B) a ~ 





B) 



A 



(6) 



FIG. 4. The top panel shows a close-up of the interesting 
region of V(4>) for the solution shown in Fig. [| The dashed 
curve show a pure V(<f>) = exp ( — X(f>) potential with A = 8.113 
for comparison. The the onset of acceleration is caused by cfi 
settling into the local minimum. The fact that the feature 
is introduced on such a small scale is due to the exponential 
factor in the potential. The lower panel shows Vp/V p and A 
(the constant curve). Where the two curves cross V' = 0. 

At one stage in this work we focused on potentials of 
the form V((j>) = exp(— X e ff<P) with the idea that X e g 
might not be absolutely constant, but could be slowly 
varying with </>. We considered forms such as X e g = A(l— 
(<j)/ B) a ) and found many interesting solutions, especially 
for moderately large values of B which make X e g slowly 
varying. For example, A = 13, B — 65, and a = 1.5 
give a solution similar to Fig ||. If this form for X e g 
were taken seriously for large <f>, then these models have 
an absolute minimum in V which cj) settles into (or at 
least approached) at the start of acceleration. But our 
expression may just represent an approximation to X e g 
over the relevant (finite) range of tf> values. Of course we 
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always can re-write Eqn. [| in terms of \ e g with X e g = 
A — ln{V p )/4>. In the end we focused on potentials in the 
form of Eqn. || because they seem more likely to connect 
with ideas from M-theory. Whatever form one considers 
for V, the concept remains the same. Simple corrections 
to pure V — exp (— A</>) can produce interesting solutions 
with all parameters O(l) in Planck units. 

We should acknowledge that we use O(l) rather loosely 
here. In the face of the sort of numbers required by 
other quintessence models or for, say a straight cosmo- 
logical constant (pa ~ 10~ 120 ) numbers like .01 and 34.8 
are 0(1). Also, the whole "quintessence" idea has sev- 
eral important open questions. Some authors argue 
that values of <f) > 1 should not be considered without a 
full quantum gravitational treatment, although currently 
most cosmologists do not to worry as long as the densities 
are << 1 (a condition our models easily meet). Another 



issue that has been emphasized by Carroll |35| is that 
even with the (standard) assumption that <p is only cou- 
pled to other matter via gravity, there still will be other 
observable consequences that will constrain quintessence 
models and require small couplings. Because in our mod- 
els <j> ~ today the tightest constraints in j35| are evaded, 
but there would still be effective dimensionless parame- 
ters w 10 required. 

Looking toward the bigger picture, a general polyno- 
mial V p will produce other features of the sort we have 
noted. Some bumps in the potential can be "rolled" over 
classically, but may produce features in the perturbation 
spectrum or other observable effects. We are investigat- 
ing a variety of cosmological scenarios with a more gen- 
eral version of V p . We are also looking at the effect of 
quantum decay processes which are relevant to local min- 
ima of the sort we consider here. We expect a range of 
possibilities depending on the nature of V p . 

In conclusion, we have exhibited a class of quintessence 
models which show realistic accelerating solutions. 
These solutions are produced with parameters in the 
quintessence potential which are O(l) in Planck units. 
Without a fundamental motivation for such a poten- 
tial, all arguments about "naturalness" and "fine tun- 
ing" are not very productive. We feel, however, that this 
work represents interesting progress at a phenomenolog- 
ical level, and might point out promising directions in 
which to search for a more fundamental picture. 
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